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We apply the non-perturbative optimized linear S expansion method to the O(N) scalar field 
| model in three-dimensions to determine the transition temperature of a dilute homogeneous Bose 

gas. Our results show that the shift of the transition temperature, AT C /T C , of the interacting model, 
compared with the ideal gas transition temperature, really behaves as -fan 1 / 3 where a is the s-wave 
scattering length and n is the number density. For N = 2 our calculations yield the value 7 = 3.059. 
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I. INTRODUCTION 



The experimental realization of the Bose-Einstein condensation in dilute atomic gases has greatly stimulated an 
enormous number of theoretical studies in this field (for recent reviews on the theory and experiments, see for instance, 
1 I Rcfs. [0,0]). Most of this interest comes from the fact that in these experiments a great deal of control can be achieved 
in almost every parameter of the system. Thus, experiments in dilute Bose gases provide a perfect ground to test 
numerous models and ideas, as for example those commonly used in quantum field theory, applied to non-relativistic 
systems. In particular, a theoretical study which has attracted some attention very recently is the determination of 
the behavior of the transition temperature in the presence of a repulsive interaction. This non trivial problem has been 
treated by different methods with different results. Taking a dilute Bose gas, with a repulsive interaction characterized 
by the scattering length parameter a, the dependence on a for the difference between the critical temperature shift 
with and without interaction (AT C /T C — (T c —T~)/T c ) is a highly controversial point, be with respect to the functional 
dependence on a or even regarding the sign. An early Hartree-Fock calculation || with a non-delta interaction gave 
a negative sign for AT C . This very same sign was also obtained by Toyoda [Q using one loop renormalization group 

and obtaining a functional dependence on a as AT C /T C = 7 (a 3 n) 1/6 , where n is the density. More recently, Huang || 
obtained the same dependence for AT C /T C , but with a positive constant 7. Griiter, Ceperley and Laloe ||, Holzmann 
and Krauth fij and Holzmann, Griiter and Laloe S] investigated the dependence of AT C numerically using Monte 

1 /3 

Carlo methods. They obtained, in the low density limit, a dependence of the type AT C /T C — 7 (a 3 n) but with 
different values for 7. More recently the Monte Carlo technique has been again applied to this problem by Prokof 'ev 
and Svistunov || and Arnold and Moore go|. These authors, who obtain 7 = 1.29 ± 0.05 and 7 = 1.32 ± 0.02, 
respectively, claim that their results are more accurate than those obtained in Refs. 



One of the reasons for the multitude of results and methods stem from the fact that at the transition temperature 
ordinary perturbation theory fails (due to infrared divergences) and we must resort to non-perturbative methods. 
Recently, there have been also used various non-perturbative methods to treat the problem of the transition temper- 
ature from an analytical way. For example, the authors in Rcf. JHJ perform a self-consistent calculation obtaining 
AT C /T C = -yan 1 / 3 , with 7 ~ 2.9. The non-perturbative 1/N method has been also used to determine the shift. Its 
leading order contribution has been evaluated by Baym, Blaizot and Zinn- Justin Q who obtain AT C /T C = -fan > , 
;_i ' with 7 ~ 2.33 for N — 2. Considering the next to leading order term, Arnold and Tomasik |p| determine a correction 
to this large-iV expansion, obtaining a value for AT C /T C which is ~ 26% smaller. 

The results obtained by Bijlsma and Stoof [Q, who used renormalization group techniques, and Baym et al. |ll]] 
were compared with the temperature transition data in the Vycor- 4 f/e system by Reppy et al. (la] . Those experiments 
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seem to give a somewhat larger value for the constant 7, as 7 ~ 5 jlfj. This value is close to 7 ~ 4.66 which is the 
one obtained in Ref. jl4| . Nevertheless some authors argue that this system would not exactly correspond to a dilute 
Bose gas of hard spheres M . 

In this paper we apply the non-perturbative linear 6 expansion (or optimized perturbation theory) |]l6| (for earlier 
references see, e.g., Refs. [[l7] ^9|) to an effective model for dilute homogeneous Bose gases. This approximation has 
been shown to be a powerful non-pcrturbative method, and sufficiently simple to use in very different applications, 
including the study of non-perturbative high temperature effects, as shown very recently in the context of finite 
temperature quantum field theory (20) as well as finite chemical potential |2l}| . The method also introduces an 
arbitrary mass parameter which prevents infrared divergence problems. Non-perturbative results are generated when 
one optimizes the theory with respect to this mass parameter. 

The paper is organized as follows. We present the method in Sec. II illustrating with an application to the pure 
anharmonic oscillator which has many similarities with the model used here to describe dilute Bose gases. The 
interpolated version of an effective model for weakly interacting homogeneous Bose gases is obtained in Sec. III. Sec. 
IV is devoted to the perturbative evaluation of density related quantities. In Sec. V we present our non-perturbative 
results for the critical temperature shift comparing with results available in the literature. The conclusions are 
presented in Sec. VI. 



II. THE LINEAR 5 EXPANSION 



The optimized linear 6 expansion (LDE) is an alternative non-perturbative approximation which has been success- 
fully used in a plethora of different problems in particle theory [^6| , p2| -[24[| , quantum mechanics 1T^J25| , |26]| , n uclear 
matter j27j, lattice field theory p^j as well as for determining the equation of state for the Ising model pj|. One 
advantage of this method is that the selection and evaluation (including renormalization) of Feynman diagrams are 
done exactly as in perturbation theory using a very simple modified propagator which depends on an arbitrary mass 
parameter. The results are optimized with respect to this parameter at the end of the calculation. The standard 
application of the LDE to a theory described by some Lagrangian density C starts with an interpolation defined by 



C s = (1 - S)C (v) +SC = C (v) + S[C - £0(77)], 



(2.1) 



where Co(rj) is the Lagrangian density of a solvable theory which can contain arbitrary mass parameters r/. The 
Lagrangian density Cs interpolates between the solvable £0(1]) (when 6 = 0) and the original C (when 6=1). To 
illustrate how the method works let us consider the anharmonic oscillator described by 
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Following the interpolation prescription given by Eq. (2.1) one may choose 
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where tt 2 = m 2 + r] 2 . The general way the method works becomes clear by looking at the Feynman rules generated 
by Cs- First, the original (j) 4 vertex has its original Feynman rule — ?6A modified to — i66X. This minor modification 
is just a reminder that one is really expanding in orders of the artificial parameter 6. Most importantly, let us look 
at the modifications implied by the addition of the arbitrary quadratic part. The original bare propagator, 



S{k) = i(k 2 -m 2 +iey 



(2.5) 



becomes 



S(k) = i(k 2 - n 2 +ie)- 1 = 



k 2 



1 - 



fc 2 - 
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(2.6) 



indicating that the term proportional to rj 2 (j) 2 contained in Cq is entering the theory in a non-perturbative way. On 
the other hand, the piece proportional to 6r] 2 (f) 2 is only being treated perturbatively as a quadratic vertex (of weight 



2 



iSr/ 2 ). Since only an i nfinite order calculation would be able to compensate for the infinite number of {—irj 2 ) insertions 
contained in Eq. (2.6), one always ends up with a 77 dependence in any quantity calculated to finite order in 5. Then, 
at the end of the calculation one sets the dummy parameter 5 to unity (the value at which the original theory is 
retrieved) and fixes 77 with the variational procedure known as the Principle of Minimal Sensitivity (PMS) |]30| which 
requires that a physical quantity P calculated perturbatively in powers of S be evaluated at the point where it is less 
sensitive to variations of the arbitrary rj. That is, one optimizes the perturbative calculation by requiring 



dP{ri) 



drj 



= . 



(2.7) 



This procedure gives fj as a function of the original parameters, including the couplings, and generates non-perturbative 
results as shown in the numerous applications cited above. 

As a warm-up for our application to the Bose gas problem we follow Belief, Garcia and Neveu j3lj evaluating the 
ground state energy density £ and the vacuum expectation value (rp 2 ) for the anharmonic oscillator. Other applications 
to this problem can be found in Refs. |T§| , |25| , p6| ] . 

By taking to = in Eq. (2.2) one obtains the Lagrangian density for the pure anharmonic oscillator (PAO) which 
cannot be treated by ordinary perturbation theory. The exact result 



g exact = A i/3 0.420804974478 . 



(2.8) 



has been obtained by Bender, Olaussen and Wang [|32| , while Banerjee et al. |33|] have obtained the exact result for 



,2^oxact = A -i/3 0.456119955748. 



(2.9) 



In quantum field theory, the ground state energy density is represented by vacuum to vacuum diagrams. The relevant 
contributions to 0{5 2 ) are |M| 
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(2.10) 



Note that the second, fourth and fifth contributions are due to the extra quadratic vertex, while all the others would 
also appear in an ordinary perturbative expansion to 0(X 2 ). Setting to = (PAO), evaluating the integrals and 
eliminating the divergent O(S ) term one obtains 
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where 
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(2.12) 



Good numerical results appear already at first order where the application of the PMS to £^{rf) yields £^ l \r]) = 
A 1/3 0.4290 at rj = (A 1.5) 1/3 . To second order the authors in Ref. Q obtain £^(77) = A 1/3 0.4210 and then carry 
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on improving this result to show convergence. The interested reader is referred to Ref. |3l| for details concerning the 
optimization procedure (selection of roots, etc). Other proofs of convergence are given in Refs. and pjjfl . 

Belief, Garcia and Neveu also investigate the vacuum expectation value (<fi 2 ) and we discuss their results here 
because this physical quantity is particularly important for our application to bosonic condensates. The perturbative 
expansion for (</> 2 ) can be obtained in different ways using standard quantum field theory methods. The authors in 
Ref. prefer to do it from the perturbative expansion for £ W recalling that 

(^)W =2 ^W. (2.13) 

Going to second order in S they optimize this quantity in two different ways. First, by applying the PMS condition 
directly to (0 2 )^ 2 ' they obtain (4> 2 )^ = A -1 / 3 0.455758. Next, they use the optimum values obtained by extremizing 
£^ 2 \ getting ((j) 2 )^ = A -1 / 3 0.454246, showing that both approaches lead to results with same order of accuracy. 

Still in the context of the anharmonic oscillator, Jones, Parkin and Winder in Ref. |34j have shown that the linear 
5 expansion applied to the calculation of dynamical evolution of (<j> 2 ), where the PMS is applied directly to this 
quantity, tracks the exact solution longer than any previous approximate methods used to study the same quantity, 
like Hartree-Fock or ordinary perturbation theory. This result also reinforces the correctness of our procedure of 
optimizing the density in this particular application to Bose condensates. 

Before applying the 8 expansion to the Bose-Einstein condensation problem let us clarify a few points regarding 
the method. Firstly, one could object to the fact that 5 is formally treated as small during the actual calculation and 
finally set to unity at the end. However, we recall that the only role attributed to this dummy parameter is to label 
the orders so that one can keep track of the extra diagrams which arise from the quadratic vertex 5r] 2 (j) 2 . 

Finally, one could ask how the LDE relates to other non-perturbative analytical methods such as the 1/N expansion. 
To see that, let us consider the same model discussed above for the case where the dynamical variables are a set of TV 
scalar fields, <\> a (a — 1, ... , TV). Proceeding as before the S expansion would give the following result for £ at O(S), 

«0) W . W J, + 4 ^«,. ( , 14) 

The application of the PMS to this quantity gives 

(TV + 2) J 1/3 



en®) 

at 



0.4290 , (2.15) 



fj= [1.5 (TV + 2)A] 1/3 . (2.16) 

Higher order contributions bring more factors of TV (more loops) making the calculation meaningless if TV is very large. 
However, this particular limit can also be properly handled by the LDE provid ed one defines g — NX declaring that 
the large-TV limit will be studied with fixed g j3j|. Using g = NX in Eq. ( [2.16 ) one sees that, in the large-TV limit, fj 



is of order N°, in terms of which Eq. (2.14) gives that £^(fj) is of order N, exactly as the leading 1/N result, as one 
can easily check. 

An important result, proved in the context of the effective potential |3(|, shows that the LDE exactly reproduces 
large-TV results in any order in 6 provided that one stays within the large-TV limit. Moreover, the LDE is sensitive 



to small-TV effects since these terms may appear in terms such as TV(TV + 2)A in Eq. (2.14). In fact, Ref. [£3[ shows 
how small- TV effects are effectively taken into account by the LDE in the context of the 1 + 1 dimensional Gross- 
Neveu model at finite temperatures where the results nicely converge, order by order, towards the exact result set by 
Landau's theorem. 

The formal relationship in between the LDE and 1/TV is investigated in detail in Refs. |23| and |36). Here, we shall 
concern ourselves with the finite TV case only. 

III. THE INTERPOLATED MODEL FOR DILUTE HOMOGENEOUS BOSE GASES 

Let us start by considering the typical model that describes a gas of interacting boson particles, described by a 
complex scalar field ip with Lagrangian density given by 
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£ = V*(x,f) + 



1 



(it 2m 



^(x,t)+^*(x, i)V(x, t) 



d 3 xV(x, i)^*(x, t)V(x - x')V'(x', *)^*(x', t) 



(3.1) 



where fi is the chemical potential. Let us take the interatomic interaction potential as being the one for a hard sphere 
gas, 



V(x-x') - — S(x-x') 

m 



(3.2) 



where a is the s-wave scattering length. 

We want to determine the deviation of the critical temperature T c , of the interacting model, in relation to the 
critical temperature for Bose-Einstein condensation for a free gas, Tq, given by the usual expression 



2n ( n 
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(3.3) 



where n is the number density of the boson gas and C(3/2) — 2.612. 

As discussed in Refs. [[[2] and Jl3|], close to the critical point we can reduce (3.1) to an effective three-dimensional 
model for the zero Matsubara frequency modes (the static modes) of the fields ip, given by the functional integration 
of the non-zero modes, obtaining an effective action defined by (/3 _1 = T) 



dr / d 3 x£ Eucl [^(x,r),^*(x,r)] -> /? / d 3 xC eS [</>(x), </>*(x)] , 



(3.4) 



where £euc1 is the Lagrangian density in Euclidean space (t = it, as usual) and with the effective action for the static 
modes, J d 3 x£ e g, being equivalent to a three-dimensional 0(2) field theory, defined by the action 



S = 



l|V^ + i^ 2 + |(^ 2 ) 2 



(3.5) 



where 4> — {4>\,4>2) is related to the original real components of tp by ip-y = (mT) 1 / 2 ^ and tp2 = (mT) 1 ' 2 ^ while r 
and u are given by 



-2m/i , u — 487ramT 



By considering the usual interpolation prescription given by Eq. (2.1) we write 

S -► S s = SS + (1 - S)S , 

where So is quadratically (exactly solvable) in the fields. 
One can choose 



where R — r + rj 1 , obtaining 
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with rj be ing an arbitrary parameter, with mass dimensions, which is fixed at a finite order in S by the PMS condition, 
Eq. (|j). Here we will optimize the physical quantity represented by (0 2 ) which, as we shall see, is directly related 
to the critical temperature shift AT C /T C . Let us first define the density for the interacting case 



n = mT((f) 2 ) u , 



(3.10) 



where, for the O(N) symmetric model, {4> 2 )u is expressed in terms of the three-dimensional dressed Green's function 
Gsip) as 
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(3.11) 



where 



G 5 (p) = [p 2 + R-6 V 2 + S 4 (p)] 1 



(3.12) 



and £5 (p) is the <fi field renormalized self energy which will be evaluated perturbatively in powers of S. 

At the critical temperature the original system must exhibit infinite correlation length, which means that at T c and 
6 = 1 (the original theory), Gj x (0) = 0. Then, one gets the relation 



-£ 4 (0) 



(3.13) 



which is just the form of the Hugcnholtz-Pincs theorem. We must stress that the choice ( |3.9[ ) respects the Hugenholtz- 
Pines Theorem at all orders in S. 



Now, by using the relation ( 3.13 ) in ( 3.11 ), one can write 
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IV. EVALUATION OF {4> 2 ) TO 0{5 2 



Expanding the above expression, in powers of 6, to 0(5) one sees that the only contribution to the self energy is a 
momentum independent tadpole diagram which is canceled by the condition on r. Then, to order <5, we obtain 
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which is u independent and cannot furnish non-perturbative results. At next order in 5 the only momentum dependent 
contribution to the self energy comes from the two loop setting sun diagram, which is of order 5 2 . Then, we obtain 
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where Y> SB (p) represents the setting sun contribution to the self energy, 
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Note that r\ acts naturally as an infrared cutoff so we do not have to worry about these type of divergences. The 
first three terms in Eq. (4.2) represent one-loop diagrams with different powers of Sr/ 2 insertions. We regularize all 
diagrams with dimensional rcgularization in arbitrary dimensions d = 3 — 2e and carry the renormalization with the 
M S scheme. So the momentum integrals are replaced by 
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where M is an arbitrary mass scale and je — 0.5772 is the Euler-Mascheroni constant. One then obtains the 0(6 2 ) 
one loop contributions 
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where we have used the expression 
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and its derivatives with respect to rj 2 to determine Eq. (4.4). The setting sun self energy diagram, with zero external 
momentum, is given by (see, for example, Ref. |37| l) 
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The momentum dependent setting sun contribution can be written as 
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We then obtain for Eq. (4.8) 
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V. THE TEMPERATURE SHIFT IN THE OPTIMIZED LINEAR 5 EXPANSION 



Using Eqs. (4.4), (4.7) and (4.10) in (4.2), we determine (4> 2 ) u at order S 2 . Note that all divergences in e cancel 
and that at order S 2 , (<j> 2 ) u is a finite quantity. One can now set S = 1 and optimize (4> 2 )u with the PMS. After that 
one sets u = in the optimized {cb 2 ) u obtaining the 5 expansion result for the critical temperature shift pH13| 
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where To is given by Eq. ([T^). At this stage it should be clear that it is preferable to optimize (cj) 2 ) 
A{(j) 2 } because the latter quantity is less rj dependent. One then obtains 
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In principle one could not sing le out one solution in favor of the other and we must be careful in choosing the 



appropriate one. Equation (5^) implies that the optimized (<fi 2 )o vanishes no matter which sign is chosen. Now (</> 2 )c 



represents the density (divided by a factor mT) in the absence of interactions, which turns out to be zero for the 
present effective theory. However, this is of no concern here since one is really interested in the difference (4> 2 ) u — (4> 2 )o, 
represented by AT C . Also, one knows that the density of the interacting gas {(4> 2 ) u ) should be smaller than that of 
the non interacting gas, which means that here one should have (4> 2 ) u < 0, which is achieved by selecting the positive 
fj. Using this in Eq. (5.1), we then get our final result, 
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Setting N = 2 in the above expression yields 



6(iV + 2)ln~ 



an3 . (5.4) 



3.059 or3 . (5.5) 



Using the 1/N expansion Baym, Blaizot and Zinn- Justin Jljj obtain AT C /T C ~ 2.33 an 1 / 3 in the leading order while 
Arnold and Tomasik |l3| obtain AT C /T C ~ 1.71 art 1 / 3 considering the next to leading order in the same approximation. 
Our result is closer to~AT c /T c ~ 2.9 an 1 / 3 obtained in Ref. |ll|] with a method which sums setting sun contributions 
in a self consistent way. These analytical results, including ours, are compared with the recent and earlier Monte 
Carlo estimates in Ref. [jlCj. Finally let us remark that the result given by Eq. ( |5.4| ) is valid only for finite N. This 
can be understood as follows. In a large N study one would have to consider uN as fixed (meaning that u — TV -1 ) 



and so, by taking N large in Eq. (4.2), one sees that the setting sun diagrams of O(N ) should be neglected since 
there are one loop diagrams of 0{N). However, these terms are linear in r\ as shown by their contribution, Eq. ( El ), 
and so the PMS does not give any meaningful result in this limit for the present model as opposed to the AO case. 
The difference arises mainly because momentum independent tadpole diagrams of O(N) are now being subtracted 
due to the Hugenholz-Pines theorem while the one loop momentum independent diagrams survive. Then, the PMS 
generates nontrivial results only by mixing diagrams which would belong to different orders in a standard large-iV 
appl ication. To have a rough idea about what is being summed one can consider the second and fourth terms in Eq. 



(4.2) to geth er with their integrated forms. Then it is clear that, apart from a numerical factor, the optimized rf given 



by Eq. (5.2) behaves as the optimized T, ss (p) — £ ss (0 ). O ne can then see that, to this order, the optimization dresses 
the simple propagator (p 2 + r/ 2 ) -1 , present in Eq. (|4.2|), with setting sun features giving non-perturbative results 



which are compatible with the ones obtained in the self consistent summation of Ref. |11[] . 

VI. CONCLUSIONS 

We conclude that our results for the critical temperature indeed reproduce the expected behavior obtained from 
other studies, which is T c ~ Tq(1 + jan 1 / 3 ). We obtain an analytical expression for the numerical coefficient 7 in terms 
of finite values of N. Our final numerical results are similar to the ones obtained with the self consistent summation 

y predicting that the numerical value of 7 is greater than the ones predicted by the 1 /N expansion at leading order 
and next to leading order |13[] . All these analytical results, including ours, have been compared with earlier and 
recent Monte Carlo results in Ref. jl0|]. It should be clear that the present calculation has been carried out to an 
order where only one two loop diagram contributes and so the quality of the approximation is hard to be inferred 
from a quantitative point of view. In fact, the purpose of the present application was just to introduce the method as 
a possible alternative to study the condensation problem. Nevertheless, one should remark that although carried out 
in a completely different fashion, our simple application seems to capture much of the features of the self consistent 
calculation performed by Baym et al. pi . 

Also, our work does not exhaust the different ways in which this method can be implemented within this particular 
problem and the possibility of further improvements is still wide open. This could be achieved by investigating 
alternative forms of implementing the method within this model, including an investigation of the best quantity to be 
optimized, and/or by pushing the calculation to higher orders. It is possible that with more refinements this method 
will generate even better numerical results with the advantage, as shown in the paper, of being considerable simpler 
and easier to use than all previous methods used to determine the behavior of T c . 

Due to its simplicity and easy implementation, we believe that the optimized 8 expansion can also be useful in other 
aspects of the theoretical study and understanding of the Bose condensation of dilute atomic gases, as determining 
the correct corrections to the energy spectrum, or in applications related to the recent studies of the dynamics of 
the Bose-Einstein condensate formation |3£| . The results of Ref. |34| are particularly motivating in the context of 
applying the optimized linear 8 expansion also to dynamical problems. Finally, we point out that Bedingham and 
Evans |39| have successfully extended the present work to the ultra-relativistic case. 
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